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C3 Abstract 

^ We consider the inverse boundary value problem in the case of discrete electrical 

Qvy networks containing nonlinear (non-ohmic) resistors. Generalizing work of Curtis, In- 

^~~^ german. Morrow, Colin de Verdiere, Gitler, and Vertigan, we characterize the circular 

,—1 planar graphs for which the inverse boundary value problem has a solution in this 

^J generalized non-linear setting. The answer is the same as in the linear setting. Our 

Q_J method of proof never requires that the resistors behave in a continuous or monotone 

r^ fashion; this allows us to recover signed conductances in many cases. We apply this to 

the problem of recovery in graphs that are not circular planar. We also use our results 
to make a frivolous knot-theoretic statement, and to slightly generalize a fact proved 
by Lam and Pylyavskyy about factorization schemes in their electrical linear group. 
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i^ 1 Introduction 

Cp Consider a network of resistors, containing several boundary nodes where current is al- 
lowed to flow in and out. If we fix the voltages at these boundary nodes, all interior volt- 

^-J ages and currents will be determined, and in particular the amount of current flowing 

px| in and out at each boundary node is determined. Thus, boundary voltages determine 

'""' boundary currents. The relationship between boundary voltages and boundary cur- 

J> rents can be encapsulated in a matrix A, the response matrix or Dirichlet-to-Neumann 

k>< matrix. 

^ Viewing the circuit as a black box whose internals are hidden, the matrix A exactly 
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describes the behavior of the circuit. It is never possible to determine the internal 
structure of a circuit from A, but if we are given the structure of the underlying graph, 
then in many cases we can determine all the resistances from A. The problem of 
recovering the resistances from boundary measurements is a discrete version of the 
inverse boundary value problem first posed by Calderon. 

We say that a graph is recoverable if the recovery problem has a solution. Working 
independently, |2j and llj completely analyzed the case of circular planar networks, 
finding a combinatorial criterion for recoverability, and finding algorithms to carry out 
the recovery when possible. In this paper, we generalize these results to allow for 
non-linear (non-ohmic) conductors. 
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In ^2]|3j we consider the forward problems, for fixed networks. Assuming that 
current at each edge is given as a continuous, monotone, zero-preserving function of 
voltage at each edge in the graph, we show that boundary voltages determine all 
internal currents uniquely. In particular, there is a well-defined Dirichlet-to-Neumann 
map from boundary voltages to boundary currents. Dually, if voltage is given as a 
continuous monotone zero-preserving function of current at each edge, then boundary 
currents more or less determine all interior currents, and there is essentially a well- 
defined Neumann-to-Dirichlet map from boundary voltages to boundary currents. 

In subsequent sections we consider the inverse problem of determining the conduc- 
tance functions at each edge from the Dirichlet-to-Neumann map or the Neumann- 
to-Dirichlet map. We are forced to assume that the conductance functions are bi- 
jective and zero-preserving, but our approach somehow doesn't require continuity or 
monotonicity, rather surprisingly. Our method of proof is apparently a variant of the 
approach used in fl], though substantial modifications are necessary, because Y-A 
transformations are no longer available. Section [9] is devoted to showing that a certain 
notion of convexity behaves well in simple pseudoline arrangements. This notion of 
convexity is used to model the propagation of information in unusual boundary- value 
problems which we use to read off conductances near the boundary, in a layer-stripping 
approach. 

Our main result is the following fact: if G is a circular planar graph for which 
recovery is possible in the usual setting with linear conductors and positive conduc- 
tances, then recovery is still possible when (sufficiently well-behaved) nonlinear and 
signed conductors are allowed. In fTTj, we apply this to two settings in which negative 



conductances naturally appear. In ^12, we restate our results from the point of view 
of Lam and Pylyavskyy's electrical linear group and consider the analogy between the 
electrical linear group and total positivity in the general linear group. 
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2 The Dirichlet-to-Neumann Map 

In this section and the next, we prove some simple facts about general networks with 
nonlinear conductors 

Definition 2.1. A graph with boundary T is a graph whose set of vertices V have 
been partitioned into a set of boundary vertices dT and interior vertices r*"*. 

We allow our graphs to contain multiple edges and self-loops, though the latter are 
completely irrelevant in what follows. If e is a directed edge in our graph, we let e be 
the edge going in the opposite direction. 

Definition 2.2. If T is a graph with boundary, a conductance network on F is a 
collection of functions 7e : R — )• M, one for each directed edge e in T, such thatjei—x) = 
— 7e(x) and such that each j^ is zero-preserving (^7e(0) = 0), monotone ('Jei^) ^ 7e(y) 
if X > y), and continuous. 

We abuse notation and use F for both the network and the graph. 

Definition 2.3. A voltage function on a graph with boundary T is a function (j) from 
the vertices ofV to M. A current function on T is a function l from the directed edges 
of F to M, such that 

L{e) = -i{e) 

and such that for every vertex v in the interior o/F, Yle ''(^) ~ ^' where the sum ranges 
over all edges leaving the vertex v. A voltage function <f) and current function l are 
compatible (with respect to some conductance network on F) if 

i{e) = -fe(.(t>{x) - 4>{y)) 

whenever e is an edge leading from vertex x to vertex y. Given such a pair, the bound- 
ary voltage function is the restriction of <j) to dV, while the boundary current function 
is the function that assigns to each vertex v E dV, the sum ^g i{e) where e ranges over 
the edges leaving from v. 



^We include these elementary results primarily because they were previously unknown within the Mathe- 
matics REU at the University of Washington (which focuses on the recovery problem in electrical networks) . 



In other words, the interpretation of the function 7e is the function that specifies 
current in terms of voltage along that edge. In the linear case, 'Jeix) = CgX, where Ce 
is the conductance along edge e. 

Given a conductance network and a set of boundary voltages, the Dirichlet bound- 
ary value problem is to find a voltage function extending the given boundary values, 
together with a compatible current function. In general, the solution to the Dirichlet 
problem may not be unique, but solutions always exist, and the current function is 
uniquely determined: 

Theorem 2.4. Let T be a conductance network and let f be any function from dT to 
M. Then there is a voltage function (p and a compatible current function l on T that 
has f as the boundary voltage function. Moreover, l is uniquely determined. 

Note that (p need not be uniquely determined. This occurs for example if every 
function 7e is the zero function and there are any interior vertices. 

Proof. We first show existence. Let R^ be the vector space of all voltage functions. 
Let W be the affine subspace of M containing only functions which agree with / on 
dT. For each directed edge e, let 
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This function is C^, nonnegative, and convex because of the properties of 7e. We also 
have (?e(0) = 0, and qe{—x) = qe{x). Moreover, if x and y are both nonnegative (or 
both nonpositive), and |x| < |y|, then 

qe{x) < qe{y). (1) 

This is easy to verify. 

For (f) G M^, we define the pseudopower Q{cl)) to be 

y^ge(0(start(e)) - (^(end(e))) 

eG-B 

where start(e) and end(e) are the start and end of e, and the sum is over all directed 
edges. Note that the terms for e and e agree, so we are essentially counting each edge 
twice. 

In the linear case, this function is the usual power, the sum over voltage times 
current for each edge. In the more general case, however, it is not this value, so we call 
it the pseudopower. 

As a sum of nonnegative convex C^ functions, pseudopower is a nonnegative convex 
C^ function on M^, and also on W. We aim to show that Q{-) attains a minimum 
on W. To see this, first let M be the maximum of |/(f)| as v ranges over dT. Let 
K be the hypercube subset of W containing only voltage functions which take values 
in [—M,M]. Then Q{-) attains a minimum on K because K is compact. Let (po be 
such a minimum on K. We claim that (po is in fact a global minimum. To see this, let 
(j) be any other voltage function in W. Let ^u : R — )■ [— M, M] be the retraction of R 



onto [-M, M] that sends (M, oo) to M and (-00, -M) to -M. Let (/>' = /i o (/>. Then 
clearly 0' still agrees with / on the boundary, so </>' G W^, and in fact (/)' G K. Moreover, 
because ^ is a monotone short map and because the qe functions satisfy Equation (II|, 
it follows easily that 

Q{<t>) > Q{ct^') > QiM- 

Since (p G W was arbitrary, (po does indeed minimize Q on all of W. 

For any directed edge e from vertex x to vertex y, let i(e) = 7e(0(x) — (/)(y)). Then l 
and (j) will be compatible current and voltage functions as long as t is a current function. 
The property that t{e) = —iie) follows easily from properties of 7e, so it remains to 
show that for every interior vertex u, 

E ^(^) = 0- 

e with u=start(e) 

This follows directly by taking the partial derivative of (5(0) with respect to (l){v), and 
using the fact that (J)q is a critical point. Specifically, the partial derivative 

^-^= E U4>{v)-<t>{^nd{e)))- E 9;(0(start(e)) -</)(!;)) 

e with ii=start(e) e with »;=end(e) 

E 7e(</'(^;) - 0(end(e))) - E 7e('/'(start(e)) - </.(f)) = 

e with u=start(e) e with i)=end(e) 

2 E ^^('^(^) - '^(end(e))) 

e with ^I=start(e) 

where the last line follows by the change of variables e <-)• e in the second sum, and the 
fact that 

— 7e((/>(start(e)) — i;^(end(e))) = 7e(</>(start(e)) — i;^(end(e))). 

This establishes existence. It remains to show that the current function i is unique. 
Suppose we had two current functions ti and t2 respectively compatible with two volt- 
age functions 4>i and 02, both of which agreed with / on dV. Define a new conductance 
network T on the same graph by letting 

7e(x) = 7e(x + 02(start(e)) - 02(end)(e)) - i2{e) 

This is a valid conductance function because 

7e (02 (start (e)) - 02(end)(e)) = i2(e) 

implies that 7e(0) = 0, and the other properties are easy to verify. One easily verifies 
that = 01 — 02 and t = ti — 12 are now compatible voltage and current functions 
for r. Moreover, the boundary values of are all zero, because 0i and 02 agreed on 
the boundary. So we are reduced to showing that if the boundary voltages vanish, all 
currents must vanish. 



Consider the directed subgraph of f in which we include a directed edge e if and 
only if t{e) > 0. Since t equals 

7e(0(start(e)) — (/>(end(e))) 

and 7e is monotone and zero-preserving, every edge e in the subgraph must also satisfy 

(/)(start(e)) > (j){end{e)) 

In particular then, our subgraph must be acyclic. We claim that it has no interior 
'sources' or 'sinks' in the following sense: there cannot be any interior node v which 
has at least one coming in and none coming out, or at least one edge going out and 
none coming in. These follow from the facts that the sum of i(e) for the edges incident 
to V is zero, and a sum of positive and nonnegative numbers cannot be zero. 

For x,y £ V, let X ^- y indicate that there is a chain of edges in the subgraph 
leading from x to y, oriented in the correct way. This relation is obviously transitive, 
and irreflexive because x — )• y implies that 4>{x) > (p{y). Suppose for the sake of 
contradiction that there is at least one edge e in the subgraph. Let e go from xq to yo. 
Then choose y with y = yo or y ^ yo, that is maximal in the sense that no y' exists 
with y ^ y'. Similarly choose a minimal x' with x' — t- xq or x' = xq. Because there 
are no internal sources or sinks, xq and yo must be boundary nodes. We also have 
Xq — > yo- Thus (p{xo) > 4'iyo), contradicting the fact that (p vanishes on the boundary. 

So this subgraph is empty. This means that for every edge e, i{e) < 0. But then 
for every edge e we also have t(e) = —t{e) < 0. So t vanishes everywhere, and ii and 
i2 agree. D 

Since the current function determines the boundary current function, it follows that 
there is a well-defined map, the Dirichlet-to- Neumann map, from boundary voltage 
functions to boundary current functions. 

3 The Neumann-to-Dirichlet Map 

This section is completely dual to the previous one. Here we assume that voltage is 
given as a suitably nice function of current, and show essentially that boundary currents 
uniquely determine boundary voltages. There are some caveats, however: a given set 
of boundary currents must add up to zero on each connected component of a graph, 
or else there are no solutions at all. Additionally, there is latitude in choosing the 
voltages: adding a constant to each value of a voltage function preserves compatibility. 
In fact, we can add a separate constant on each connected component of T. 

Definition 3.1. IfV is a graph with boundary, a resistance network on T is a function 
that assigns to each directed edge e a function p^ : R ^ R. which is 

• zero-preserving: pe{0) = 

• monotone: Pe{x) > pe{y) if x'>y 

• continuous 



such that pe{—x) = —pe{x). 

Definition 3.2. A voltage function (/) and current function t on T are compatible 
(with respect to some resistance network on F) if 

Pe{i{e)) = (t>{x) - 4>{y) 

whenever e is an edge leading from vertex x to vertex y. Given such a pair, the bound- 
ary voltage function is the restriction of (p to dT, while the boundary current function 
is the function that assigns to each vertex v E dV, the sum ^^ i{e) where e ranges over 
the edges leaving from v. 

In other words, the interpretation of the function pe is the function that specifies 
voltage in terms of current along that edge. In the linear case, Pe{x) = r^x, where r^ 
is the resistance (reciprocal of conductance) along edge e. 

The dual to Theorem |2.4| is the following: 



Theorem 3.3. Let f he a real-valued function on dV, summing to zero on each con- 
nected component ofT. Then there is a voltage function (f) and a compatible current 
function l with f as the boundary current function. Moreover, given two such solutions 
(pi and (j)2, (pi ~ 4'2 is constant on each connected component ofT. 

Proof. We first show existence. Let W be the affine space of current functions on 
r which have / as their boundary currents. This space is nonempty because of the 
requirement that / sum to zero on each connected component - we leave the proof of 
this as an exercise to the reader. For each edge e, let 



qe{x) = / Pe{t)dt 

Jo 

As before, these functions are nonnegative, convex, and C^ and satisfy Equation ([I|. 
For t eW, let Q{i) be the sum 

e 

Then as before, Q(-) is convex, nonnegative, and C^. 

Let M = X^^ggr 1/(^)1 ^^"^ ^^^ ^ be the set of all current functions l £ W for which 
there are no directed cycles of edges ei,e2, . . . , e„ with t(ej) > for all i. We claim 
that K is compact, because \i{e)\ < M for i £ K and e any edge. This is true because 
the acyclicity of l implies that we can partition the vertices of F into a set of vertices 
'upstream' from e and a set of vertices 'downstream' from e, and the total current 
flowing through e must be at most the total current flowing into the upstream set, 
which is at most M. So K is bounded, and clearly closed. Thus Q attains a minimum 
on K, at some to- 

As before, we wish to show that to is a global minimum. To see this, take any 
i G W \ K. Then l contains a directed cycle of edges along which current is flowing 
in a positive direction. Letting m be the minimum of these positive current values, 
create another current function l' by subtracting m from the current along each of these 
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edges (and adding m to the currents along the edges going in the opposite directions). 
Then clearly t' is still a current function, and it has the same boundary currents as 
L. Moreover, since the current did not change in sign along any edge, and did not 
increase in magnitude along any edge, it follows by the analogue of Equation M that 
Q{i') < Q{i^)- Now l' may have more cycles, and so may not yet be in K. However, 
by repeating this procedure at most finitely many times, we get an l" G K with 
Q{i) > Q{i") > Q{l,o). The process must terminate after finitely many steps because 
the number of edges along which no current is fiowing increases at each step. Since i 
was arbitrary, lq is indeed a global minimum. 

Take i = lq. Then we need to produce a voltage function (j) for which 

4>{x) - 4>{y) = Pe{^{e)) 

whenever e is an edge from x to y. This is possible exactly when 

5]pe,We.)) = (2) 

i 

for every directed cycle ei, 62, . . . , Cn- To see that Equation Q holds, take our cycle 
ei, 62, . . . , e^, and let ip be the current function on T which takes the value 1 along 
each of the Cj, — 1 along each ei, and elsewhere. Then t/) is a current function with 
boundary current zero, so L + ttp G W for every i G M. Then the fact that i is a global 
minimum of Q{-) on W implies that 

dQJL + tjj) 
dt 
vanishes at t = 0. But by an easy direct calculation, this gives Equation ([2]). 

Next we show that (p is essentially uniquely determined. Let <j)i and 02 be two 
solutions, respectively compatible with ti and L2- As in the previous section, we can 
shift the resistance functions and produce a new network which has <j)i — 02 as a voltage 
function compatible with the current function ii — 12- So without loss of generality, 
/ is identically zero, and we need to show that if and i are compatible voltage and 
current functions with boundary current zero, then is constant on each connected 
component. 

Similar to the previous theorem, we take a subgraph containing all edges e for which 

0(start(e)) > (/)(end(e)) 

Since voltage is given as a monotone function of current, this implies that i(e) > 
too. As before, there can then be no cycles and no internal sources and sinks. But in 
fact, since the boundary current is also zero, there can be no sources or sinks on the 
boundary either. But any finite directed acyclic graph without sources or sinks must 
be empty. Therefore, for every edge e we must have 

0(start(e)) < (/)(end(e)) 
But then applying this to e, we see that 

0(start(e)) > 0(start(e)) 
too. Thus (j) must be constant on each connected component. D 



Let A be the vector space of real- valued functions / on dT with the property that 
for every connected component C of F, 

vecndT 

Then it is easy to see that there is a well-defined Neuniann-to-Dirichlet map A ^ A 
from boundary voltage functions in A to boundary current functions in A. In the case 
of a conductance or resistance network where the conductance or resistance functions 



are bijections, it follows by combining Theorems 2.4 and 3.3 that we have a well-defined 
bijection between boundary voltage functions in A and boundary current functions in 
A. In other words, the Dirichlet-to-Neumann map and the Neumann-to-Dirichlet map 
are essentially inverses of each other. Also, boundary voltages uniquely determine 
interior voltages in this case, so the solution to the Dirichlet problem is unique. 

4 The Inverse Problem 

Given the Dirichlet-to-Neumann map or the Neumann-to-Dirichlet map, can we deter- 
mine the resistance or conductance functions? Without putting additional restrictions 
on our conductance or resistance functions, the answer is no in almost all cases, be- 
cause of the following problem. Consider the Y-shaped network of Figure [T] If we put 
conductance functions of arctanx on two of the three edges, then the total current 
through the third edge can never exceed vr. Consequently, if we put a conductance 
function of 73(3;) = x on the third edge, this yields the same Dirichlet-to-Neumann 
map as the alternative conductance function 



73(2;) 



X if — TT < X < TT 

■n if X > TT 
^ — vr if X < — vr 



Similar arguments work as long as any interior vertices exist. In the dual case where 
voltage is given as a function of current, a similar problem occurs in any graph that 
contains cycles. 

To resolve this difficulty, we instead consider networks where the conductance func- 
tions or resistance functions are bijections: 

Definition 4.1. A bijective network is a graph F with boundary, and a function that 
assigns to each directed edge e inV, a conductance function 7e : M — )• M, which satisfies 
the following properties: 

• 7e(0) = 

• 7e is a bijection for every e 

• le{-x) = -le{x) 

We say that the network if monotone if every 7e is order-preserving (x > y ^ 7e{-c) > 
7e(y)J, and we say that the network is linear if every 7e is linear (of the form 7e(x) = 
CeX for some Cg € Mj. 
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Figure 1: A graph for which the general recovery problem is impossible. The dark middle 
vertex is an interior vertex, while the light vertices are boundary vertices. 

For a fixed graph T with boundary, a bijective conductance structure on T is a 
choice of the functions je- 

The definition of compatible current and voltage functions and boundary values is 
the same as for conductance networks. 

Note that we do not require the 7e to be continuous. In the case of monotone bijec- 
tive networks, however, the je are all continuous, so by combining Theorems |2.4| and 
Theorems 13.31 we have well-defined Dirichlet-to- Neumann and Neumann-to-Dirichlet 
maps. In the general case of bijective networks, current may not be uniquely deter- 
mined by voltage, and vice versa, so we use a relation instead of a function: 

Definition 4.2. If T is a bijective network with n boundary nodes, then the Dirichlet- 
Neumann relation is the relation A C M" x M" consisting of all pairs {u,f), where u 
is the boundary voltage of a voltage function 0, / is the boundary current of a current 
function i, and 4> and t are compatible. 

The inverse problem is the problem of determining the conductance functions of a 
bijective network given the Dirichlet-Neumann relation and the underlying graph. 

Definition 4.3. A graph T with boundary is strongly recoverable if no two distinct 
bijective network structures on T have the same Dirichlet-Neumann relation. T is 
weakly recoverable if no two distinct monotone linear bijective network structures on 
r have the same Dirichlet-Neumann relation. 

In other words, a graph is strongly recoverable if the inverse problem has a solution, 
in the sense that the conductance functions are uniquely determined by the Dirichlet- 
Neumann relation. A graph is weakly recoverable if this holds with the assumption 
that the conductance functions are monotone linear. Weak recoverability is the type 
of recoverability considered in [2] and [1] . 

Definition 4.4. Let T be a graph with boundary. We say that T is circular planar if it 
can be embedded in a disk D with all boundary vertices on dD and all interior vertices 
in the interior of D. 
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Figure 2: Covoltages on the left, and the corresponding currents on the right. Arrows point 
in the direction of current flow. 



Equivalently, T is circular planar if the graph obtained by adding an extra vertex 
V and connecting v to every vertex in dV is planar. 

Generalizing results of p] and [I], we will determine exactly which circular planar 
graphs are recoverable, showing that for circular planar graphs, weak recover ability 
and strong recover ability are equivalent. The techniques of [2j used linear algebra and 
are no longer applicable in a non-linear setting, while the techniques of [l] will require 
modification because the Y-A transformation is no longer available. 



5 Covoltages 



Working with voltages and currents is somewhat asymmetric because voltage functions 
live on vertices and are unconstrained, while current functions live on edges and are 
constrained by Kirchhoff 's Current Law. To make matters more symmetric, it will help 
to work with covoltages rather than currents. 

Definition 5.1. If T is a circular planar graph with boundary (with a fixed planar 
embedding), then a covoltage function on T is an arbitrary real-valued function on the 
faces ofV. If ip is a covoltage function, the associated current function is the current 
function which assigns a value of ilj{fi) — V'(/2) to a directed edge with face /i on its 
left and /2 on its right. 

See Figure [2] for an example. We summarize the key facts about covoltage functions 
in the following lemma: 

Lemma 5.2. 

• The associated current function of a covoltage function is a current function (it 
satisfies KCL). 

• If ip is a covoltage function, with associated current function l, then the boundary 
current of i at a boundary vertex v is ip{fi) — ip{f2), where /i and /2 are the faces 
on the left and right sides of v (from v facing into the graph). See Figurel^ 
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Figure 3: A boundary vertex v and the two faces /i and /2 on either side. 

• A covoltage function is determined by its associated current function up to a 
constant of integration. In other words, if two covoltage functions have the same 
associated current function, then their difference is constant. 

• Every current function is associated to a one covoltage function. 

Only the last of these is nontrivial. It follows from well-known facts about planar 
graph duality, and we leave its proof as an exercise to the reader. 

From this lemma we see that covoltage functions and current functions are essen- 
tially equivalent. We will generally identify them, saying that a voltage function (j) is 
compatible with a covoltage function ip if (p is compatible with the associated current 
function of ip. 

More explicitly, a voltage function </> is compatible with a covoltage function ijj if 
at every edge e from vertex vi to V2, with face /i on the left and /2 on the right, 

4>{vi)-Hv2) = leWfl)-i^{f2))- 

From the lemma, we also see that boundary covoltages determine boundary cur- 
rents, and conversely boundary currents determine boundary covoltages up to an ad- 
ditive constant. Consequently, the Neumann-Dirichlet relation contains the same in- 
formation as the voltage- covoltage relation H which can be defined by 

^ = {(<A|<9V) V'I^-f) : is a voltage function and ^ is a compatible covoltage function} 

where dV and dF are the sets of vertices and faces of F along the boundary. 

The relation between H and A can be described explicitly as follows: let vi, f i,V2, f 2, ■■■ ,Vn, fn 
be the boundary vertices and faces of F, in counterclockwise order around the boundary 
of F. Then we have 

(5, /i) G H ^^ {g{vi),g{v2), ..., g{vn), /i(/„)-/i(/i), Hfi)-h{f2), ..., /i(/„-i)-/i(/n)) G A, 

and conversely, (xi, X2, • . . , x„, yi, 2/2, . . . , y„) G A if and only if yi + 7/2 H h y„ = 

and {g, /i) G H where g{vi) = Xi and h is any function on the set of boundary faces with 
h{fn) — ^(/i) = 2/1) h{fi) — h{f2) = 2/2) and so on. Such a function h exists because 
2/1 + 2/2 + • • • + yn = 0, and the choice of h doesn't matter because H has the property 
that {g, /i.) G H if and only if {g, h + k) £ Xi. 

Because the voltage-covoltage relation contains the same information as the Neumann- 
Dirichlet relation, we can replace the Neumann-Dirichlet relation with the voltage- 
covoltage relation in the definition of weak and strong recoverability 
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Figure 4: A sample medial graph. This one is rather degenerate; one of the geodesies closes 
back up on itself. 

6 Medial Graphs 

A graph's recoverability can be determined by considering an associated medial graph: 

Definition 6.1. A medial graph is a piecewise- smooth Jordan curve C C M^^ called the 
boundary curve and a collection of piecewise- smooth curves (^ geodesies j in the inside of 
C , each of which is either closed, or has both endpoints on C , satisfying the following 
properties: 

• Every intersection of a geodesic with C is transversal. 

• Every intersection between two geodesies or between a geodesic and itself is transver- 
sal. 

• There are no triple intersection points of geodesies or the boundary curve. 

The geodesies partition the inside of the boundary curve into connected components 
cahed cells. Those along the boundary curve are called the boundary cells. The interior 
vertices of a medial graph are all the points of intersection between geodesies (including 
self- intersections) . The boundary vertices are the intersections of the geodesies with 
the boundary curve. 

We draw the boundary curve with a dashed line. The medial graph of Figure |4] has 
three geodesies, eleven cells, eight interior vertices, and four boundary vertices. 

To any circular planar graph F (with a fixed embedding in a disk D), we can 
construct a medial graph as follows: place a vertex along each edge of F, and connect 
two vertices if they are on consecutive edges around one of the faces of F. Also place 
two vertices along each arc of dD between consecutive boundary nodes of F, and 
connect these to the vertices on the closest edges. Figure [5] shows this whole process. 
The resulting vertices and edges can then be decomposed as a union of geodesies in a 
unique way. The result is a medial graph with at least \dT\ geodesies, and boundary 
curve dD. 

Definition 6.2. A coloring of a medial graph is a function from the cells to the set 
{black, white} such that directly adjacent cells have opposite colors. Here we say that 
two cells are directly adjacent if they share an edge in common. 
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Figure 5: The process of converting a circular planar graph (top left) into a medial graph 
(bottom left). The natural coloring of this medial graph is shown on the bottom right. 
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Figure 6: Medial graphs that cannot come from circular planar graphs (regardless of which 
way they are colored). 
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Figure 7: The four medial-graph cells surrounding an edge of the original graph. 

Every medial graph has exactly two colorings. A medial graph that comes from 
a circular planar graph F has a canonical coloring, with the cells that come from 
vertices colored black, and the cells that come from faces colored white. The process 
of converting a circular planar graph into a colored medial graph is almost reversible. 
In particular, a circular planar graph T is determined (up to isotopy) by its associated 
colored medial graph. However, some degenerate medial graphs, like those in Figure |6| 
do not come from circular planar graphs. 

Colored medial graphs therefore generalize circular planar graphs, and we can gen- 
eralize the inverse boundary value problems to colored medial graphs in the following 
way. If M is the medial graph of F, then the black cells of M correspond to the ver- 
tices of F, while the white cells of M correspond to faces of F. So a voltage function is 
equivalent to a function on the black cells of M, while a covoltage function is equivalent 
to a function on the white cells of M. The compatibility condition then says that if w, 
X, y, z are four cells that meet at an interior vertex corresponding to a directed edge 
e, then 

ip{z) - 4^{x) = -YeiHw) - 0(y)) 

if w, X, y, z are in counterclockwise order around e, and w is the cell at the start of e, 
as in Figure [7j From the medial graph point of view, we might as well combine ^ and 
(j) into one function. 

Remark 6.3. To avoid a definitional nightmare, we will henceforth pretend that con- 
ductance functions 7e are always odd. Equivalently, we pretend that 7e = 7e- None 
of the following results will use this assumption, but they are much harder to state 
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correctly without it! For example, in the next definition, we would need to define some- 
thing like an "oriented interior vertex" which is tricky to define correctly in the case 
where geodesies loop back on themselves. 

Definition 6.4. Let M be a colored medial graph. A conductance structure on M 
assigns to every interior vertex v an odd bijection 7^, : M — )■ M. A conductance structure 
is monotone if every 7^ is order-preserving, and linear if every 7t, is linear. 

When M comes from a circular planar graph F, a conductance structure on M is 
equivalent to a bijective conductance structure on T, and the conditions of monotonicity 
and linearity are equivalent for M and for F. 

Definition 6.5. Let M be a colored medial graph and ^ be a conductance structure on 
M. A labelling of (M, 7) is a function (p from the cells of M to M such that for every 
interior vertex v in M , 

ct>{z) - (Pix) = -f,{ct>{w) - cl>{y)) (3) 

where w, x, y, z are the cells in counterclockwise order around v, with w and y colored 
black. The boundary data of a labelling (p is the restriction (P\q of 4> to the boundary 
cells. The boundary relation H^ of 7 is the set 

■^7 = Wld '■ (p is a labelling of (M,^)} 

of all possible boundary data of labellings of (M, 7) . A colored medial graph M is 
strongly recoverable if the function 7 1— )■ H^ is injective, and M is weakly recoverable 
if it is injective when restricted to monotone linear 7. 

Note that strong recoverability implies weak recoverability. In the case where M 
comes from a circular planar graph, it is not difficult to see that H^ contains the same 
information as the voltage-covoltage relation of (F,7), and in particular, F is weakly 
(strongly) recoverable if and only if M is. 

If M' is the colored medial graph obtained by reversing the colors of M, then M 
is weakly (resp. strongly) recoverable if and only if M' is weakly (resp. strongly) 
recoverable. In particular, recoverability of a medial graph doesn't depend on the 
choice of the coloring. We won't use this fact in what follows, so we leave the proof to 
the reader. 

Definition 6.6. A medial graph is semicritical if the following hold: 

• No geodesic intersects itself 

• // gi and g2 are geodesies, then gi and 52 intersect at most once. 

If in addition, every geodesic starts and ends on the boundary curve (rather than being 
a closed loop), we say that the medial graph is critical. 

Critical medial graphs are essentially equivalent to simple pseudoline arrangements. 
The main result we are working towards is the following: 

Theorem 6.7. Let M be a medial graph. Then the following are equivalent: 
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• M is strongly recoverable. 

• M is weakly recoverable. 

• M is semicritical. 

In particular, this gives a criterion for determining whether a circular planar graph is 
weakly or strongly recoverable. 

The equivalence of weak recoverability and semicriticality for medial graphs coming 
from circular planar graphs was shown in [2] and [I] by arguments involving linear 
algebra and determinantal identities. It therefore comes as a surprise that this remains 
true in the highly nonlinear cases considered here. It is also interesting that for circular 
planar graphs, weak and strong recoverability are equivalent. It is unclear whether this 
remains true for non-planar graphs. 

7 Lenses, Boundary Triangles, and Motions 

In this section, we review facts about medial graphs from Is], and show why a weakly 
recoverable medial graph is semicritical. Almost all of the proofs of this section are 
based on those in |^. We use unusual terminology, however. 

Definition 7.1. A boundary digon or boundary triangle is a simply- connected bound- 
ary cell with two or three sides, respectively. If g is a geodesic which does not intersect 
the boundary curve or any other geodesies and which is a simple closed curve, then we 
call the set of cells inside g a circle. // there is just one cell, we call it an empty circle. 

Definition 7.2. If M is a critical medial graph and g is a geodesic, then since M is 
critical, g divides M into two connected components, by the Jordan curve theorem. We 
call such a connected component a side. If g and h are two geodesies, then gUh divides 
M into four pieces, which we call wedges. The boundary of any wedge consists of a 
segment of the boundary curve, which we call the base, and two segments of g and h, 
which we call the legs. We call the intersection of g and h the apex of a wedge. We 
say that a wedge W is a semiboundary wedge if at least one of the two legs contains 
no interior vertices other than the apex. A boundary zone is one of the following: 

• A semiboundary wedge 

• The side of a geodesic which intersects no other geodesies. 

Theorem 7.3. Let M be a critical medial graph. If M has more than one cell, then 
M has at least one boundary digon, or at least three boundary triangles. 

Proof. We first prove the following: 

Lemma 7.4. Every boundary zone contains a boundary triangle or a boundary digon. 

Proof. Let Z he a boundary zone. Without loss of generality, Z contains no smaller 
boundary zones. I claim that no geodesic crosses the boundary of Z. This is trivial if 
Z is the side of a geodesic which intersects no other geodesies. In the other case, Z is a 
semiboundary wedge. Say Z has apex a and legs ab and ac. Without loss of generality, 
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no geodesies eross ab, and it remains to show that no geodesies eross ac. Suppose a 
geodesic g crosses ac at an interior vertex d. Without loss of generahty, d is as close to 
c as possible, and in particular, there are no interior vertices between c and d. After g 
enters Z at d, it must do one of three things: 

• Exit Z by passing through ab. 

• Exit Z by passing through ac. 

• End on the boundary curve. 

The first of these is impossible because there are no interior vertices on ab. The second 
is impossible because it would mean that g crosses ac in two places. The third is 
impossible because it would mean that there was a smaller semiboundary wedge inside 
Z, as in Figure [Sj contradicting the choice of Z. 

So no geodesic enters or exits Z. I claim that there are in fact no geodesies in 
Z. Let g be any geodesic in Z. If g intersects no other geodesies, then one side of 
5 is a smaller boundary zone inside Z, a contradiction. Otherwise, there is at least 
one interior vertex along g. Let e be one which is as close to the boundary curve 
as possible, as in Figure [9} Then e is the apex of at least two semiboundary wedges 
contained inside Z, a contradiction. 

Therefore, Z contains no geodesies in its interior. It follows that Z must be a single 
cell. Therefore Z is either a boundary triangle or a boundary digon. D 

We next, show the existence of at least one boundary triangle or boundary digon. 
By the lemma, it suffices to find at least one boundary zone. Since M has more than 
one cell, it has more than zero geodesies. Let g he a geodesic. If g intersects no other 
geodesies, then both sides of g are boundary zones, so we have at least one boundary 
digon or boundary triangle. Otherwise, let t; be a boundary vertex of g which is 
as close as possible to one of the endpoints of 5. If /i crosses g at v, then half the 
wedges formed by g and h are semiboundary wedges, so at least one boundary digon 
or boundary triangle exists. 

If a boundary digon exists, then we are done. Otherwise, we have a boundary 
triangle. This boundary triangle is in fact a wedge formed by two geodesies g and 
h. Two of the other wedges formed by g and h are semiboundary wedges, so this 
gives two more boundary triangles or boundary digons. This completes the proof of 
Theorem I7.3[ D 



Definition 7.5. A medial graph is circleless if it has no empty circles. 

Theorem 7.6. A semicritical circleless graph M is critical. 

Proof. Given the definition of critical and semicritical, it suffices to show that if M is 
circleless, then it has no circles. Let C be a circle which minimizes the number of cells 
inside C. Since all the intersections of geodesies are transversal, if g is any geodesic, 
then g intersects (the boundary of) C an even number of times. This can only happen 
in a semicritical graph if g intersects (the boundary of) C exactly zero times. It follows 
that every geodesic in M is either inside C or outside C. But every geodesic inside C 
must itself be a circle, contradicting the choice of C. Therefore C is an empty circle, 
and M is not circleless. D 
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Figure 8: If there are any interior vertices along leg ac, then abc is not a minimal semibound- 
ary wedge. 




Figure 9: If any interior vertex e exists along g, then e is the apex of at least two semibound- 
ary wedges inside abc. 
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Figure 10: An empty 2-lens, an empty 1-lens, and a nonempty sharp 2-lens. 

Definition 7.7. If g is a geodesic, we can think of g as either a continuous function 
on [0, 1] or a continuous function on M with period 1. If [a, 6] is a subinterval of the 
domain of g, and g is infective on both half-open intervals (a, 6] and [a, 6), we call the 
curve g{[a,h]) a geodesic segment with endpoints g{a) and g{h). 

Definition 7.8. Let v be an interior vertex. A 1-lens with pole at v is the interior 
of a Jordan curve that is a geodesic segment with endpoints v and v. Let v and w be 
interior vertices. A 2-lens with poles at v and w is the interior of a Jordan curve that 
is a union of two geodesic segments, each having endpoints at v and w. For n = 1,2, 
we say that an n-lens is sharp if for every pole v, only one of the four cells that meet 
at V is inside the n-lens. We say that an n-lens is empty if it contains only one cell. 

Because every medial graph can be colored, no cell can be on both sides of a 
geodesic. Prom this it is easy to see that empty lenses must be sharp lenses. See 



Figure 10 for examples of lenses. 



Theorem 7.9. If M is a medial graph which is not semicritical, then M must contain 
a sharp 1-lens or sharp 2-lens. 

Proof. Since M is not semicritical, we either have a geodesic which self-intersects, or 
two geodesies which intersect too many times. First suppose that some geodesic g 
intersects itself. Then we have g{s) = g{t) for some s ^ t, and in the case where g is 
a closed curve and we think of 5 as a function on M/Z, we want s — b ^ Z. Choose a 
pair {s, t} for which all this holds, minimizing |s — 1|. Without loss of generality s < t. 
I claim that g{[s, t]) is a geodesic segment. In the periodic case, clearly |s — t| < 1, or 
else we could replace s with s -|- 1 and make |s — i| smaller. If g fails to be injective on 
{s,t], then there exist s < ri < r2 < t with g{ri) = g{r2) and |ri — r2| < |s — t| (and 
ri — r2 ^ Z in the periodic case), contradicting the choice of s,t. A similar argument 
shows that g is injective on [s,t). So g{[s,t]) is a geodesic segment with endpoints at 
g{s) and g{t). In particular, it is a 1-lens, with pole at f = g{s) = g{t). 

If this 1-lens fails to be sharp, then some geodesic h at v enters into the 1-lens. It 
eventually exits the 1-lens. If it does so without crossing itself, then we have a smaller 



2-lens or 1-lens, as in Figure 11 Otherwise, we can repeat all of the above argument on 
h (or rather, on h restricted to the part between v and where h first exits the 1-lens), 
and get a smaller 1-lens, and continue inductively, until we get a sharp 1-lens. 

Next, suppose that no geodesic intersects itself. Then M can only fail to be semi- 
critical if two geodesies g and h intersect more than once. So there exist s,t in the 
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Figure 11: A non-sharp 1-lens must contain a smaller 2-lens or 1-lens. A third possibility, 
not pictured here, is that the geodesic crosses itself before entering and exiting the original 



1-lens. Then one repeats the argument of Theorem 7.9 inductively. 




Figure 12: If sharpness fails at pole f , then geodesic g has no means of exiting the 2-lens. 
(The pole at w needn't be sharp in this arrangement.) 



domain of h with h{s) and h{t) in the range of g, and s ^ t. As before, we can take 
s — f ^ Z in the periodic case. Choose such an {s,t} with \s — t\ minimized. Then as 
before, |s — 1| < 1 in the periodic case. By choice of {s, t}, h{r) cannot lie on g for any 
r between s and t. It follows immediately that the geodesic segments between h{s) and 
h{t) along h and g form a 2-lens, with poles at t> = h[s) and w = h{t). To see that this 
2-lens is sharp, suppose that more than one of the cells around v lies inside the 2-lens. 
Since h ^ g and the intersection at v is transversal and the boundary of the 2-lens is 
a Jordan curve, exactly three of the cells lie inside, and we have the configuration of 



Figure 12, Then g continues on from v into the interior of the 2-lens. It has no way of 
escaping, however: it cannot intersect itself (by assumption), and it cannot intersect 
h{\s,t\) (by choice of {s,i}). This is a contradiction. D 



Definition 7.10. A Y-A motion is the transformation shown in Figure \Tl 

If our medial graphs come from circular planar networks, this is nothing but a 



standard Y-A transformation as in Figure 14 



Theorem 7.11. If M is a medial graph which is not semicritical, then we can apply 
some combination of the following operations to reach a medial graph which contains 
an empty 1-lens or 2-lens: 
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Figure 13: A Y-A motion in a medial graph. 





Figure 14: A Y-A transformation. 



• A Y-A motion 

• Removal of an empty circle. 

Proof. By the previous theorem, there is at least one sharp lens. We prove by induction 
on n that if M is a medial graph containing a sharp lens which contains n cells, then 
M can be converted into a medial graph which contains an empty lens. The base case 
where n = 1 is trivial, because then M already has an empty lens. Suppose n > 1, 
and we have a sharp lens L with n cells. If there are any smaller sharp lenses inside L, 
we are done by induction. Suppose that L contains no smaller sharp lenses. Since the 
boundary of L is a Jordan curve, we can make a new medial graph M' from the inside 
of L, with the boundary of L as its boundary curve. I claim that M' is semicritical. 
Otherwise, it contains a sharp lens, by the previous theorem, contradicting the fact 
that L contains no smaller sharp lenses. If M' contains an empty circle, then we can 
remove it, and continue by induction. So we can suppose that M' is circleless. Then by 



Theorem 7.6, M' is critical. Therefore, it contains at least three boundary triangles, 
or at least one boundary digon. If there are at least three boundary triangles, then at 
least one of them does not touch a pole of L. Therefore it corresponds to an actual 
triangle in M, and we can carry out a Y-A motion, pushing the intersection outside 
of L, while preserving the lens L. This decreases the number of cells in L, so we can 
continue by induction. Finally, suppose that there is a boundary digon. If it does 
not touch one of the poles of L, it corresponds to an empty lens, and we are done. 
Otherwise, it corresponds to a triangle in M, and we can push it outside of L by a 
Y-A motion, and then continue by induction. D 

In the monotone linear case, a conductance structure on a medial graph is merely 
a function that assigns a conductance c^ > to each interior vertex v of M. The 
conductance functions are of the form jvix) = CyX. 
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Figure 15: A Y-A motion in a colored medial graph. 

Lemma 7.12. The following transformations do not effect the boundary relation of 
a medial graph. We assume that the vertices involved in the transformations have 
monotone linear conductance functions, hut the rest of the medial graph need not. 

1. Removing a boundary digon. 

2. Removing an empty circle. 



3. A motion as in Figure 15. Here, the new conductances should he given by 

ClC2 



Cl2 



Cl3 



C23 



Ci + C2 + C3 
C1C3 

Ci + C2 + C3 
C2C3 

Ci + C2 + C3 



in one direction, and 



c\ 



c-2 



C3 



C12C13 + C12C23 + C13C23 

C23 
C12C13 + C12C23 + C13C23 

Cl3 

C12C13 + C12C23 + C13C23 



C12 



in the other direction. 
4. Removing an immediate self-loop. 



5. Replacing an empty lens with a single intersection, as in Figure 16. In the top 



row of Figure 16, the conductance of the new vertex should he 

C1C2 



Cl +C2 



while in the bottom row it should he 



Cl +C2 



Proof. These all follow by direct calculation. The third rule is just the standard Y-A 
transformation, while the fifth consists of the standard rules for resistors in series or 
parallel. D 
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Figure 16: Replacing an empty lens with a single intersection. 



The Y-A transformation is reversible, while the process of replacing a series or 
parallel configuration with a single resistor throw away information. This immediately 
yields the following conclusion: 

Corollary 7.13. 

• Adding or removing a boundary digon does not affect weak or strong recoverahility. 

• Adding or removing an empty circle does not affect weak or strong recoverahility. 

• A Y-/S. motion does not affect weak recoverahility. 

• Any medial graph containing an immediate self-loop or an empty lens is not weakly 
or strongly recoverable. 



This yields half of Theorem 6.7 



Corollary 7.14. If M is a medial graph and M is not semicritical, then M is not 
weakly recoverable. 



Proof. Combine Corollary |7. 13 and Theorem 7.11 



D 



For the other direction, we will use the following property of boundary triangles. 
Theorem 7.15. Let M he a medial graph with a boundary triangle abc, and let M' be 



the medial graph obtained by uncrossing the boundary triangle, as in Figure 17. Let a 
be the apex of abc. If we fix the conductance function 'ja, then the boundary relation of 
M' is determined by the boundary relation of M . 



Proof. Let w,x,y,z be the four cells around a, as in Figure 17 By ([s]), if (f> is any 
labelling of M', then (p extends to a labelling of M uniquely by setting 

4>{w) = (piy) ± jaiH^) - 4>{y)) 

where the it depends on the coloring. It follows that if H and r! are the boundary 
relations of M and M' , then 



(...,p,g,r, ...) G H' 
so H and H' determine each other. 



(...,p,g±7a(p-r),r, ...) e H, 



D 
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Figure 17: "Uncrossing" a boundary triangle. 



In light of this, we can reduce Theorem 6.7 to one condition: 



Theorem 7.16. Suppose the following statement is true: if M is any critical colored 
medial graph, and abc is a boundary triangle of M , with apex a, then 7^ is determined 
by the boundary relation of M . Then Theorem\6.7\ is true. 



Proof. We have already seen that weak recoverability implies semicriticality. Clearly, 
strong recoverability implies weak recoverability. It remains to show that semicriticality 
implies strong recoverability. Let M be a semicritical medial graph. Suppose M is not 
strongly recoverable, and take M to have as few cells as possible. By Theorem |7. 13 



we 



can remove empty circles from M without affecting strong recoverability. Therefore, 



M is circss, and by Theorem 7.6, it is critical. Similarly, M cannot have a boundary 



digon, by Theorem |7.13[ If M has no geodesies, then M is strongly recoverable because 
the inverse boundary problem requires the recovery of no information. Otherwise, M 



must have a boundary triangle by Theorem 7.3, Let M' be the medial graph obtained 
by uncrossing the boundary triangle. By assumption, we can find the conductance 
function at the apex of the boundary triangle, and by the previous theorem, this in 
turn yields the boundary relation of M' . Uncrossing a boundary triangle does not break 
criticality, so M' is critical, and by choice of M it is therefore strongly recoverable. So 
given the boundary relation of M' , we can find all the remaining conductance functions 
of M. Therefore M is strongly recoverable, a contradiction. D 

It remains to show how to recover boundary the conductance functions of boundary 
triangles. We will do this by setting up special boundary value problems which allow us 
to read off the conductance function at a boundary trianglejj The main difficulty will 
be showing that our boundary value problem is well-posed, having a unique solution. 

8 The Propagation of Information 

To accomplish this, we will need a way of modelling the propagation of information 
through a medial graph. We begin by defining partial labellings of medial graphs. 

Definition 8.1. If M is a fixed medial graph, a cellset is a set of cells in M. If M is 
the medial graph associated to a bijective network T, and S is a cellset, then a labelling 
of S is a function </> : S — )• M which satisfies ^ whenever w, x,y, z S S. 



^In 1 2 , conductances of boundary triangles are recovered using linear algebra. Since we allow nonlinear 
conductances, this approach won't work here. 
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Figure 18: The wavy cellset has rank six, because it contains seven cells, and there is one 
interior vertex completely surrounded by S (the circled one). 

Note that the notion of a labelling depends on the conductance functions. 

Next we turn to combinatorial definitions that depend solely on M, not even on its 
coloring. These will be our model for the propagation of information. For motivation, 
see Theorem 18/71 below. 



Definition 8.2. Let S be a cellset. Then define rank(S') to be \S\ — v{S) where \S\ is 
the number of cells in S, and v{S) is the number of interior vertices for which all four 
surrounding cells are in S. 

The rank of S estimates the number of degrees of freedom of a labelling of S, since it 
counts the number of cells (the number of variables) minus the number of constraining 
equations. See Figure [18] for an example. 

Definition 8.3. If S is a cellset and x is a cell not in S, then we say that Su{x} is a 
simple extension of S if there is an interior vertex v of M for which three of the cells 
around v are in S and x is the fourth. If there does not exist another interior vertex 
v' with the same property, then we say that S U {x} is a safe simple extension. 



Lemma 8.4. Let S be a cellset and x be a cell not in S. Then rank(S'U{x}) < rank(S') 
if S U {x} is a simple extension of S, and rank(5 U {x}) = rank(S') + 1 otherwise. 
Moreover, rank(S'U{x}) = rank(5') if and only if SL){x} is a safe simple extension of 
S. 



Proof. Obviously \S U {x}\ 
need to show that 



|5| + 1. So letting v{T) = rank(S') — |5| as above, we 
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Figure 19: The top right is a safe simple extension of the top left, while the bottom right is 
a (non-safe) simple extension of the bottom left. 
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Figure 20: The forbidden configuration in a closed set. A cellset is closed exactly when this 
configuration never occurs. 

• If S" U {x} is not a simple extension of S, then v{S U {x}) = v{S) 

• If S" U {x} is a safe simple extension of S, then v{S U {x}) = v{S) + 1. 

• If S" U {x} is a simple extension of S that is not safe, then v{S U {x}) > v{S) + 1. 

But these are all obvious from the definitions of simple extensions and safe simple 
extensions. D 

Definition 8.5. A cellset S is closed if it has the property that for every interior vertex 
V of M , if three of the four cells around v are in S, then so is the fourth. The smallest 
closed set containing S is the closure of S, denoted S. 

The closure S exists because the intersection of closed cellsets is closed. 

Lemma 8.6. A cellset S is closed if and only if it has no simple extensions. If S' is 
a simple extension of S, then S' = S. Also, there is a chain of cellsets 

S = SiCS2C---CSn = S 

where S'j+i is a simple extension of Si for each i. Moreover rank(S') < rank(S'), with 
equality if and only if Si^i is a safe simple extension of Si for each i. 

Proof. The first claim is obvious. The second claim follows because S <Z S' implies 
that S C S', but conversely S' is clearly a subset of S and so 5' C S* because S' is the 
smallest closed subset containing S. 

For the third claim, inductively define Sj+i to be any simple extension of Si, stop- 
ping at Sn which has no simple extension. Then Sn is closed by the first claim, and 
Sn = Sn = S hy the second claim. 

For the final claim, we know by the previous lemma that rank(S'j+i) < rank(S'j) 
for every i. Then by transitivity rank(S') < rank(S'), and if equality holds then 
rank(6'j+i) = rank(S'i) for every i, so that each step in the chain is a safe simple 
extension. D 

The significance of rank and closure lies in the following: 



Theorem 8.7. Let S be a cellset and (j) he a labelling of S. If S' is a simple extension 
of S, then (j) extends to a labelling of S' in at most one way, and if S' is a safe simple 
extension of S, then (p extends in exactly one way. Consequently, cp extends to a 
labelling of S in at most one way, and i/rank(5) = rank(S'), then it extends in exactly 
one way. 

Proof. First note that Equation ([s]) above has the property that any three of 4>{w) , (f>{x) , (^{y) , 
and ^(z) uniquely determine the fourth, because we stipulated that 7e be bijective. 

If S' = SL){x}, then by definition of simple extension, there is some interior vertex 
V of the medial graph such that x is one of the four cells around v, and the other three 
are already in S. Then (j){x) is completely determined by ([3]) at the edge corresponding 
to V. Moreover, if S'U {x} is a safe simple extension, then v is the unique vertex which 
determines 4>{x), so (p{x) is uniquely determined. Thus a labelling of S extends to a 
labelling of 5' in at most one way, and at least one way if S' is a safe simple extension 
of 5. 

The statements about extending (p to S then follow easily from the previous lemma. 

D 

Definition 8.8. A cellset S is safe i/rank(S') = rank(5). 

We have just shown that if S is safe, then a labelling of S extends to a labelling of 
5 in a unique way. 

9 Convexity in Medial Graphs 

From now on we will assume that our medial graph M is critical, so that each geodesic 
starts and ends on the boundary, no geodesic intersects itself, and no two geodesies 
intersect more than once. The results of this section are completely independent of 
conductances, and are entirely combinatorial in nature. Since critical medial graphs 
are equivalent to simple pseudoline arrangements, our results can all be viewed as 
statements about pseudoline arrangements, but we will prefer to work with medial 
graphs (which contain a boundary curve) to simplify some proofs. 

If 5^ is a geodesic in a medial graph, the Jordan curve theorem implies that g divides 
the medial graph into two parts, which we call pseudo halfplanes. We can identify these 
with cellsets. 

Definition 9.1. A cellset S is convex if it is an intersection of zero or more pseudo 
halfplanes. The convex closure of a cellset S is the smallest convex set containing S, 
i.e., the intersection of all half-planes containing S. 

Since half-planes are closed, so are all convex sets. We will roughly show the 
converse. In doing this, we will show that our notion of convexity has many of the 
usual properties of convexity. Many of our results would be obvious if our medial 
graphs were always equivalent to medial graphs drawn with straight lines, but this is 
not always the cascQ 



On page 101 of 11 , Ringel gives a simple pseudoline arrangement which cannot be drawn with straight 



lines. It is obtained by perturbing the Pappus configuration. 
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Definition 9.2. // a and h are two cells in a medial graph, we say that a and b 
are adjacent if they share an edge. A path between a and b is a sequence of cells 
a = xq, xi, . . . ,Xn = b for n > 0, with Xi adjacent to Xj+i for each i. The number n is 
the length of the path. A minimal path between two cells is a path of minimal length. 
The distance between a and b, denoted d(a, b), is defined to be the length of a minimal 
path between a and b. A cellset S is connected if for every a and b in S, there is a 
path between a and b containing only cells in S. The connected components of S are 
the maximal connected subsets of S. 

Each cellset is the disjoint union of its connected components. The medial graph 
as a whole is always connected, so that d(a, b) is well-defined. 

Lemma 9.3. If S is a connected cellset, then so is S. 



Proof. Because S can be obtained from 5 by a series of simple extensions (by Lemma 8.6 ) 
it suffices to show that a simple extension of a connected cellset is connected, which is 
obvious. n 

Definition 9.4. Let g be a geodesic, and a, b be cells. We say that g separates a and 
b if a and b are on opposite sides of g. We let d'(a,6) denote the number of geodesies 
which separate a and b. If S is a cellset, we say that g separates cells of S if there exist 
a,b £ S on opposite sides of g. 

Clearly, d'{a,b) <d(a,6). 

Theorem 9.5. Let S be a nonempty convex cellset. Then there is a critical medial 
graph M' whose cells are the same as the cells in S. In particular, S is connected and 
the boundary of the closure of the union of the cells in S is a Jordan curve. Additionally, 
the geodesies in M' are in one to one correspondence with the geodesies of M which 
separate cells of S. 

Proof. If i/ is a pseudo halfplane, then we can make a medial graph out of the cells 
in H, as shown in Figure [2T| This new medial graph will still be critical. Moreover, if 
H' is any other pseudo halfplane, then H' H is a pseudo halfplane in the new medial 
graph, by criticality. We show by induction on r that if M is a critical medial graph 
and S is the intersection of r pseudo halfplanes in M, then S is the set of cells of a 
critical medial graph. The base case where r = is trivial. Otherwise, write S as 
H n S', where H is a half-plane, and S' = 01=1 ^i ^^ pseudo halfplanes H^. Then we 
can consider H as a medial graph M' , and H f] Hi is a pseudo halfplane of M' for each 
i. By induction, S = H Ci S' = 01=1 i^ ^ ^i) is the set of cells of a critical medial 
graph, and we are done. We leave the second claim as an exercise to the reader. D 

Corollary 9.6. If a and b are two cells of a critical medial graph M, then d'(a,6) = 
d(a,6). If c is a third cell, then c is on a minimal path between a and b if and only if 
c is in the convex closure of{a,b}. 
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Figure 21: A pseudo halfplane (left), considered as a medial graph in its own right (right). 

Proof. We proceed by induction on d'(a, b). If a and b are the same cell, then clearly 
d(a,6) = d'(a,b) = 0. Assume therefore that a ^ b. Let S be the convex closure of 
{a,b}. By the theorem, S is connected. Since S is connected, there is a path from o 
to b in S, and in particular there is some cell a' £ S which is adjacent to a. Let g be 
the geodesic which separates a' from a. Then g also separates a from b. Otherwise, 
the pseudo halfplane H on the side of g containing a and b would not contain a', so 
a' ^ S. Since g separates a from b and o from a', it follows that d'(o', b) = d'(a, 6) — 1. 
By induction, d'(a, 6) = d'(a', 6) + 1 = d(a', 6) + 1. Therefore, 

d'(a, 6) < d(a, 6) < d(a, a') + d(a', 6) = d(a', 6) + 1 = d'(a, 6), 

so the desired equality holds. 

For the second claim, note that c is on a minimal path between a and b if and only 
if d{a,b) = d(a,c) + d(c, 6). Define ag{x,y) for cells x,y and geodesic 5 to be 1 if 5 
separates a from b, and otherwise. Then the first claim of this corollary implies that 
d{x,y) = ^„o-g(x,y). Now any geodesic which separates a from b must separate o 
from c or c from b, so 

cTgia, b) < ag{a, c) + ag{c, b). 

Summing over (7, we see that d(a,6) < d(a,c) + d(c, 6), with equality if and only if 
(7g(a, b) = (Tg{a, c) + o"g(c, b) for every geodesic g. Thus c is on a minimal path from a 
to 6 if and only if for every geodesic g, either g does not separate any of a, 6, c, or g( 
separates a from 6. Equivalently, c is on a minimal path from a to 6 if and only if c is 
in every half plane which contains both a and c. But this is just the convex closure of 
{a, 6}. D 

Corollary 9.7. Every cell in a critical medial graph has a boundary that is a Jordan 
curve. 

Proof. Let a be a cell, and let S be the convex closure of {a}. If b is any adjacent cell 
to a, then some geodesic g separates a from b. Then there is a pseudo halfplane H 
which contains b but not a, so that b ^ S. Thus S contains no cells adjacent to a. But 
by Theorem 9.5, S is connected, so S* = {a}. It then follows that a itself must have the 



shape of a medial graph, and in particular its boundary must be a piecewise smooth 
Jordan curve. D 
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Figure 22: The desired setup in the proof of Theorem |9.9[ Here, we can enumerate the 
neighbors of a in order. In particular, they occur contiguously. 



Definition 9.8. Let (xi, X2, . . . , x„) and {yi,y2, ■ ■ ■ , Vn) be two paths from a to b. We 
say that (yi, . . . , y„) is a simple deformation of (xi, . . . , x„) if there is some 1 < i < n 
such that Xj = yj for j ^ i, and Xj_i, Xj+i,Xj and yi are the four cells surrounding 
Sonne internal vertex. We say that (2/1,2/2, • • • , yn) is a deformation of (xi, . . . , Xn) if it 
is obtained by a series of zero or more simple deformations. 

Theorem 9.9. Ifa^ b are cells, then all minimal paths between a and b are deformations 
of each other. 



Proof. We proceed by induction on d(a,6). The base cases where d(a, 6) 
obvious, so assume that d(a, 6) > 1. 



or 1 are 



By Corollary 9.6, the union of all the minimal paths between a and b are the convex 



closure of {a,b}, which by Theorem 9.5 is a medial graph itself. So we can assume 



without loss of generality that the convex closure of {a, b} is everything, and so every 
geodesic separates a from b. 

Let S be the set of all cells immediately adjacent to a. We claim that we can 
enumerate the elements of S* as si, ,52, . . . , s„ in such a way that for each i, the three 
cells a. Si, and Si-\-i meet at a vertex (with a fourth cell). This setup is shown in 



Figure 22 This is obvious if a is not a boundary cell (which probably never happens) . 
If a is a boundary cell, then this is still clear, unless a meets the boundary of the 
medial graph in two or more disjoint boundary arcs, as in Figure [23} But in this case, 
we can add a new geodesic connecting two of these boundary arcs, and throw away the 
side of the geodesic that is opposite g, making a smaller medial graph. Any geodesic 
discarded by this operation could not have separated a from b, and could not have been 
present. Consequently, the discarded part could not have contained any additional cells 
or geodesies, implying that the two boundary arcs we connected were in fact the same, 
a contradiction. 

So we can number the neighbors of a as si, S2, . . . , s„. Any minimal path from a 
to b must begin by moving from a to one of the Si. By induction, we know that all 
minimal paths from Sj to b are equivalent through deformations. So it suffices to show 
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Figure 23: A cell whose neighbors are not contiguous. This only occurs if the cell touches 
the boundary more than once, as a does here. 



for each i that at least one path from a to 6 through Si is a deformation of a path from 
a to 6 through Sj+i. 



To see this, let c be the fourth cell meeting a, Si, and Sj+i, as in Figure 24 If we 
let Qi and gi+i be the geodesies separating a from Sj and Sj+i, respectively, then both 
Qi and Qi+i separate a from b. Consequently, d(c, b) = d(a, b) — 2. Thus if we take a 
minimal path from c to b, and add on (a, Sj) or (a, Sj+i) to the beginning, we will get 
minimal paths from a to b. But these two alternatives are deformations of each other, 
so we are done. D 

The following theorem gives several conditions equivalent to convexity. 

Theorem 9.10. Let S be a cellset. The following are equivalent: 

(a) S is closed and connected. 

(b) S is connected, and for every a,b £ S, if one minimal path between a and b is in 

S, then every minimal path between a and b is in S. 

(c) If a,b £ S, then every minimal path from a to b is in S. 

(d) If a,b are adjacent cells separated by geodesic g, and a £ S but b ^ S, then every 

cell in S is on the same side of g as a. 

(e) S is convex 
Proof. 



(a)=^(b) By Theorem 9.9, it suffices to show that if 5 is closed, and P is a path in S, 
then every deformation of P is in S. This is trivial, by definition of closed. 
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Figure 24: The setup of the rest of the proof of Theorem 9.9 
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(b)=^(c) Let S satisfy (b). We first show that S is closed. Let w,x,y,z be four cells 
in S that meet around an interior vertex, in counterclockwise order, Suppose that 
x,y,z £ S. If x = z, then the cell x contradicts Corollary |9.7[ So x / z. And since 
medial graphs are always two-colorable, x and z are not adjacent. Consequently 
x — )• tf — )• 2; and X ^ y ^ z are both minimal paths from x to z. Since the former 
is in S, the latter must also be in S, showing that S is closed. 
Using this, we show that for every a,b £ S, there is a minimal path between a 
and b in S, which suffices to show (c) given (b). Since S is connected, between 
any two a,b G S there is at least one path. Call a path S -minimal if it has 
minimal length among all paths that lie in S. Then it remains to show that every 
S'-minimal path is truly minimal, since an S'-minimal path always exists between 
two cells in S. 

We proceed by induction on the length n of our S'-minimal path. The base cases 
where n < 2 are trivial. So suppose that xq —>• xi —>••••—)■ a;„ is an S'-minimal 
path for n > 2. Then the path xi — t- • • • — t- x„ is also S-minimal, so by induction 
it is truly minimal. In particular, 

d(xi,x„) = n- 1. 

Now let g be the geodesic which separates xq from xi. If x„ is on the same side 



of g as xi, then by Corollary 9.6, d(xo, x„) = d(xi, Xn) + 1 = n and we are done. 
So we can assume that x„ is on the same side of g as xq. 

Similarly, the path xq —)•••• —)• x„_i is minimal. Since it begins with a step 
across the geodesic g, Xn-i must be on the opposite side of g from xq. This 
makes x„_i and x„ on opposite sides of g, so they must be separated hy g. In 
other words, the path xq — ?• • • • — t- x,i begins and ends by stepping across g. Now 
the path xi —)•••• —)• x„_i is S-minimal, so by induction it is truly minimal. But 
then property (b) implies that every minimal path between xi and x^-i is in S. 
In particular, the path from xi to x„_i that proceeds directly along the side of 



geodesic g, as in Figure 25, must be in S. This path is minimal because it steps 
over no geodesic twice, by criticality. 

Consequently, every cell on the far side of g from xi to x^-i is in S, as are the 
cells Xq and x„ which are opposite xi and x„_i, respectively. Thus we have the 



setup of Figure 25 But then the fact that S is closed easily implies that every 
cell on the near side of g from xq to x„ is in S. In particular, there is a truly 
minimal path from xq to x„ in S, so every S-minimal path from xq to x„ must 
be truly minimal. 

(c)=^(d) Let S satisfy (c), and suppose that a and b are adjacent cells, with a £ S and 
b ^ S. Let H be the pseudo-halfplane containing a but not b. We wish to show 
that S (^ H. Suppose for the sake of contradiction that c G S and c ^ H. Then by 



property (c) and Corollary 9.6, the convex closure of {a, c} is in S. We claim that 
b G {a, c} C S, a contradiction. Otherwise, there must be some pseudo-halfplane 
H' containing a and c but not b. But then the geodesic which defines H' must 
separate a from 6, so H' must be H or its complement. Either way, H' cannot 
contain both a and c, since a £ H and c ^ H. 
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Figure 25: The path in question must begin and end by crossing g. By induction, the path 
that proceeds directly from xi to Xn~i along the side of g must be in S. Once all these cells 
are in 5, the cells on the near side of g must also be in 5, because S is closed. This shows 
that a true minimal path from xq to Xn is in S. 



(d)= 



>(e) Let S satisfy (d). If S is empty, the result is trivial, so suppose that S is 
nonempty. Call g a defining geodesic of S if there are adjacent cells a, h separated 
by g, with a ^ S and b ^ S. In such a circumstance, call the half-plane on the 
side of g containing a a defining halfplane. Then condition (d) implies that S is 
contained in every defining halfplane. Let S' be the intersection of all defining 
halfplanes. We have 5" C S", and S' is convex. It remains to show that S = S' . 



(e) = 



Otherwise there exists some cell in S' \ S. By Theorem 9.5 5' is connected, so 
there is some path from a cell in S to a cell in S' \ S. In particular, there exists 
some pair of adjacent cells a,b £ S' with a £ S and b ^ S. But then the halfplane 
containing a but not & is a defining halfplane of S , so b ^ S' , a contradiction. 

>-(a) Let S be convex. Then Theorem 9.5 shows that S is connected. Convex sets 



are always closed, because half-planes are closed and intersections of closed sets 
are closed. 

D 

Corollary 9.11. Let S be a cellset. Then S is closed if and only if each connected 
component of S is convex. 

Proof. It is easy to check that S is closed if and only if each connected component of S 
is closed. The result then follows by the equivalence (a) <^=^ (e) of Theorem 9.10 D 



Next we return to the notion of rank: 

Theorem 9.12. If S is a convex cellset, then rank(S) is one plus the number of 
geodesies g which separate cells of S. 



Proof. By Theorem 9.5 it suffices to show this in the case where S is the entire medial 
graph. Then we need to show that the rank of S is one plus the number of geodesies. 
Let c count the number of cells, Vi count the number of interior vertices, Vb count the 
number of boundary vertices, and e count the number of edges including the segments 
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of the boundary curve. Interior vertices have degree 4 while boundary vertices have 
degree 3, so that 

2e = 4vi + 3ffe. 



Since all cells in a critical medial graph are homeomorphic to disks (by Corollary 9.7), 
Euler's formula applies, yielding 

3 

c + Vi + Vb = e + l = 2vi + -Vb + 1. 

Thus 

rank(5) = c - Vi = -Vb - ft, + 1 = — + 1. 

But Vb/2 counts the number of geodesies, because each geodesic has two endpoints. D 

Recall that a cellset is "safe" if its closure has the same rank as it. The main 
statement that we will need to set up our boundary value problems is the following: 

Lemma 9.13. Let S be a safe cellset with connected closure. Let b ^ S be a cell 
adjacent to a cell in S. Then S L)b is a safe cellset with connected closure. 



Proof. Note that S'U{5} = S L) {b} by general abstract pr ope rties of closure opera- 



tions. Since b is adjacent to S and S is connected. Lemma 9.3 shows that S U {6} 
S U {b} is connected. 

Since both S and S U | b| are closed and connected, they are convex by Theo- 
rem 



9.10 



So Theorem 9.12 tells us that the ranks of S and S U {b} are determined by 



the number of geodesies which cut across them. Now any geodesic which cuts across 
S will also cut across S* U {6}, so 



rank(S) < rank(5 U {6}). 



But in fact at least one geodesic which cuts across S U {6} does not cut across S, 
specifically the geodesic g which separates b from its neighbor in S. This geodesic 
cannot cut across S by condition (d) in Theorem 



9.10 



So we have 



rank(5) + 1 < rank(5 U {b}). (4) 

Now 

rank(S') = rank(S') 

because S is safe. And since b ^ S, S D {b} is not a simple extension of S, so that 

iank{S U {b}) = iank{S) + 1 
by Lemma |8.4[ And by Lemma |8.6[ 



rank(5 U {b}) < rank(5 U {b}) = rank(5) + 1 = rank(5) + 1. 
Combining this with Q, we see that equality must hold, so in particular 



rank(5 U {b}) = rank(5 U {b}), 
so 5" U {6} is safe. D 

We will use this lemma to inductively build up safe sets of boundary cells. 
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S\T 




Figure 26: The setup of Lemma 10.1 Given boundary cell b and adjacent geodesic g, we 
can find sets S and T of boundary cells having the desired properties. 



10 Nonlinear Recovery 



To complete the proof of Theorem |6.7[ we need to show how to recover the conductance 
functions of boundary triangles in critical medial graphs. 

Lemma 10.1. Let M be a critical medial graph, b be a boundary cell, and g be one of 



the geodesies by b, as in Figure 26. Then there are two sets of boundary cells T, S such 
that 

• S is a safe cellset: rank(S) = rank(S'). 

• rc5. 

• beS\T. 

• T is the halfplane on the far side of g from b. 

• S is the entire medial graph. 

Proof. Let H be the half-plane on the far side of g from b. Number the boundary cells 
in H ai,02, . . . , a„ in consecutive order. Let Tq = {ai}, and then recursively define 
Tj+i to be Ti U {ojil where ji is the minimal j such that aj ^ Tj. Stop once Tj contains 
every Uj . This gives us a chain of sets of boundary cells 



Ti c Ts c • • • c r„ 



By choice of Oj. , we must ha ve a,. 



-1 G Tj for each i. Consequently, Oj^ is adjacent to 
a cell in Tj. Then by Lemma 9.13, we see that if Tj is safe with connected closure, so 
is Tj+i = Tj U {aj.}. Since Ti = {a} clearly has Ti = Ti, Ti is safe with connected 
closure. Thus every Ti is safe with connected closure. 

Let T = Tn- Then T is safe with connected closure, and T contains every boundary 
cell in H. Since T<ZH,T<^H = H. But conversely, if any half-plane H' contains T, 
then its defining geodesic g' must start and end on the far side of H' , so by criticality 
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it cannot cross (7, and H C H' . Consequently the convex closure of T contains H. But 
by Theorem 9.10, T is already convex, so we see that H QT. So -^= T. 

TnU{b}. Now b iT 



9.13 



again, l^+i is 



We then let Tn+i = T„ U {6}. Now b ^ T = H, so hy Lemma 
still safe with connected closure. Continuing on inductively we can keep adding more 
boundary cells until we have a set S which is safe, has connected closure, contains 
TU {6}, and has every boundary cell in its closure. But S is convex by Theorem 9.10 
and any convex set containing every boundary cell must be the whole medial graph, 
since no halfplane contains every boundary cell. 

D 



See Figure 27 for an illustration of all these steps in a small example. 



Using this we prove the other direction of Theorem 6.7 



Proof (of Theorem 6.1). By Theorem 7.16, we only need to show that conductance 



functions of boundary triangles are recoverable. Let 6 be a boundary triangle, and let 



g be one of the two geodesies which define b. Find cellsets T and S as in Lemma 10.1 



Fix boundary values of x at cell 6, and on S'\ {6}. Since S is safe and has closure the 



entire medial graph, it follows by Theorem 8.7 that these boundary values uniquely 
determine a labelling cf) of the entire medial graph. The restriction of </) to T is a labelling 
of r, and so is the constant zero function. Since these extend the same labelling of T, 
Theorem 8.7 shows that (j) must identically vanish on T. By Lemma 10.1, T contains 



a neighbor a of 6, and the cell c that diagonally touches a. Letting d be the fourth cell 
that meets at the apex of 6, we have from Equation ([3]) 

<l^{d) = c^id) - m = ±7^(./.(a) - </.(c)) = ±7^(x) 

where the signs depend on the orientation of e, the order of a, 6, c, d around their 
meeting point, and whether or not e is a boundary spike. But we can determine (t){d) 
from the boundary relation H, since d is a boundary cell. It follows that we can read 
off ib7^(x) directly from H. D 



11 Applications of Negative conductances 



Even though our original motivation for Theorem 6.7 was the recovery of networks 
that had nonlinear monotone conductance functions, monotonicity was never needed 
in the proof. In fact, we could have allowed our currents and voltages to take values 
in an arbitrary abelian group. This result is rather strong, and one wonders what sort 
of applications it might have. 

One interesting implication is that we can allow linear conductance functions 7(3;) = 
ex arbitrary c G C \ {0}, and recovery will be possible. While these nonpositive 
conductances might seem arbitrary, there are a couple of circumstances where they 
come up naturally. 
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Figure 27: A small example of the proof of Lemma |10.1[ By a series of safe simple extensions, 
we first fill up the desired pseudo halfplane (the one below the darkened geodesic), and then 
proceed to fill up the rest of the medial graph. Arrows show the propagation of information. 
Here, g is the darkened geodesic and the boundary cell h is the one above the left end of g. 
S consists of all four boundary cells below g, while T additionally contains h and the next 
two boundary cells after h in clockwise order. 
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Figure 28: A complete graph on four vertices, and a critical circular planar graph. Both have 
six edges, and generically, they are interchangeable, after an appropriate birational change 
of variables. 



11.1 Removing a mild failure of circular planarity 



One circumstance is a generalization of the Y-A transformation of Section|7| | Consider 



the circular planar network shown on the right side of Figure 28 , and suppose that all 
conductors are linear. By a direct calculation, one can show that this network is 
electrically equivalent to the network on the complete graph i^4 on the left side of 



Figure 28 with the conductances specified as follows: 



\12 



a + 



cd 







a 


Al3 = 


b + 


ce 
a 


Al4 = 


cl 
a 




A23 = 


de 
a 




A24 = 


dl 
a 




A34 = 


^ 





a 



where a = c + d + e + f, provided that c + d + e + f. Conversely, if the Ajj are positive 
numbers, then the 1^4 on the left side of Figure [28] is electrically equivalent to the graph 



12 



in 1995 and in Michael Goff's REU 



^The following idea appears in Konrad Schr0der's REU paper 
paper [5] in 2003. The latter paper claims to have a proof that signed conductances can be recovered in 
critical circular planar networks, provided that the Dirichlet-Ncuniann relation is a function. However, the 
proof seems flawed to me. 
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Figure 29: A K4 with positive conductances can be turned into one of the four critical 
circular planar graphs of the top row, possibly introducing negative conductances in the 
process. The four associated medial graphs are shown on the bottom row. 



on the right side of Figure 28 with conductances given as follows: 



a = X 



A 



c 
d 



12 



13 



A14A23 

A34 
A14A23 



24 



qX 



14 



A24A34 

q 

A34 
g 

A24 



A 



where 



23 



A14A23 + A23A24 + A23A34 + A24A, 



^24^34 



So we have a way of transforming any K^ (with positive conductances) into an 



electrically equivalent network on the graph on the right side of Figure 28, possibly 
introducing nonpositive conductances in the process. Now from the equations for 
a, . . . , /, we see that only a or 6 can be zero or negative. A conductance of zero is 
equivalent to a deleted edge, so any K4^ (with positive conductances) can be turned 
into one of the four graphs of Figure 29 (with conductances in M \ {0}). These graphs' 



associated medial graphs are also shown in Figure 29 



Like a Y-A transform, we can apply this transformation to a i^4 embedded locally 



within a larger graph. For example, the graph of Figure 30 (with positive conduc- 



tances) will be electrically equivalent to one of the graphs of Figure p^ (with nonzero 
conductances). Now these resultant graphs are circular planar, and it is easy to check 



that they are critical. Consequently, the original graph of Figure 30 is (weakly) recov- 
erable - as long as the response matrix provides enough information to distinguish the 
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Figure 30: If all the conductances in this graph are positive, it will be electrically equivalent 



to one of the four graphs in Figure 31 , all of which are recoverable. 



four graphs of Figure 31 This is implied by the following lemma, which generalizes 
results in [s]. 

Definition 11.1. Let Mi and M2 be two critical medial graphs, embedded on the same 
disks with boundary vertices in the same locations. Then we say that Mi and M2 are 
equivalent if for every boundary vertex v, the geodesies qmi of Mi and (7^2 of M2 
beginning at v end at the same boundary vertex. 

Lemma 11.2. Let Mi and M2 be two colored critical medial graphs with conductances, 
with the same boundary vertices. If Mi and M2 have the same boundary relation, then 
Ml and M2 are equivalent. 

Proof. Let w be one of the common boundary vertices of Mi and M2. Let gi be the 
geodesic in Mi with one endpoint at w, and let Vi be the other endpoint of gi. Suppose 
for the sake of contradiction that vi ^ V2. Let H be the half-plane in Mi on the side of 



gi containing V2, and let b be the cell of Mi that is by w but not in H, as in Figure 32 
By Lemma 10.1 we can find sets of boundary cells T C S, with b £ S \T, S and T 
safe, T = H, and S equal to the whole medial graph Mi. So we can set boundary 
values of on T, 1 at b, and whatever we like on 5" \ {b}, and there will be at least one 
labelling of Mi with these boundary values. In fact, since T = H, the labelling must 
take the value zero on all of H, so we can even assert boundary values of zero on every 
boundary cell in H. 

In particular, we see that for Mi, there is at least one set of valid boundary values 
that takes the value 1 at b, and at every boundary cell in H. We claim that this 
fails in M2, a contradiction. To see this, let Q be the set of boundary cells in M2 
corresponding to the boundary cells of H. That is, Q consists of all the boundary cells 
in M2 along the boundary arc from w to vi that passes through V2. Let b' be the cell 
of M2 corresponding to b. If Mi and M2 had the same boundary relation, then there 
should be a labelling of M2 which takes the value at every cell in Q, and 1 at b' . 
But this is impossible because b' £ Q, so that boundary values of zero on Q should 
force a value of zero at b' . To see that b' € Q, note first that Q is connected, so Q 
is connected and thus convex. But b' is directly across 52 from a cell of Q, so if any 
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Figure 31: These four graphs are recoverable. Furthermore, they are electrically distinguish- 



able from each other, by Lemma 11.2 



44 




V2 




Q 



(a) Ml (b) M2 

Figure 32: Two nonequivalent medial graphs. 





Figure 33: Turning a "star" into a complete graph. 



geodesic separates b' from Q, it would be §2, which is impossible because cells on both 
sides of (72 are in Q, specifically the two cells around V2- So no geodesic separates b' 
from Q, and b' G Q. 

So the assumption that vi ^ V2 was false. Consequently gi and (72 have the same 
start and endpoints. D 



By Lemma 11.2 and an examination of the medial graphs of Figure 31 , we see that 



all are electrically distinct, and consequently the original graph of Figure 30 is weakly 
recoverable. 



In contrast, the graph of Figure 34a is not recoverable because in one of the cases, 



it is equivalent to Figure 34b, which is not recoverable. To verify this, note that we 



can turn a "four-star" into a K4 as in Figure 33, without producing negative or zero 
conductances. This map is invertible (because the "four-star" is recoverable), so if 
Figure 34a were recoverable, then Figure |34b would also be recoverable, but it is not 
since its medial graph clearly contains a lens. 

However, we can say that Figure 34a] is generically recoverable, because generically 
it is equivalent to Figure 34c which is weakly recoverable. In fact, the map from 



conductances on Figure 34a to response matrix is birational. 
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(b) 




Figure 34: Of these three graphs, only |34c| is recoverable, though |34a| is generically recover- 
able. 



11.2 Knot Theory 

Another place where negative conductances naturally appear is in knot-theoretic stud- 
ies of tangles. As explained in [6], we can associate a medial graph (with conductances) 
to a picture of a knot-theoretic tangle, by viewing the strands in the tangle as the 
geodesies, and assigning a conductance of ±1 at each crossing, with sign determined 
by which of the two strands is on top. The electrical behavior of the network is invari- 
ant under the Reidemeister moves, making it a topological invariant. In the case of 
knots, which can be thought of as tangles with no loose ends, the resulting electrical 
network has no boundary vertices, so no useful information is obtained this way. In the 
case of tangles with four loose ends, the resulting electrical network has two boundary 
vertices, so its electrical behavior is summarized by a single number, the resistance. 
For the well-known family of rational tangles (see [7]), the resistance is a complete in- 
variant of the tangle. In particular, up to equivalence there is a unique rational tangle 
for every element of Q U {oo}. 

Combining facts from p\ and u\ with our results on signed conductances, one easily 
obtains the following novelty theorem: 

Theorem 11.3. Fix some pseudoline arrangement (i.e., critical medial graph) with n 
pseudolines. View the geodesies as strands in a knot diagram or tangle diagram. At 
each crossing in the original pseudoline arrangement, substitute a rational tangle not 
equivalent to or oo. Interpret the resulting diagram as a tangle with 2n loose ends. 
Then given the initial pseudoline arrangement and the resulting tangle, one can say 
which rational tangle was substituted at each vertex. 



As an example, the two tangles in Figure |35[a) are both rational tangles not equiv- 
alent to or oo. Substituting them into the pseudoline arrangement of Figure 35 ^b) 
in two different ways, Theorem 11. 3| implies that the two tangles in Figure 36 
equivalent. 



are not 



It might be interesting to consider whether Theorem 11.3 remains true if "rational 
tangle" is replaced by some broader class of 2-strand tangles. 
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(a) (b) 

Figure 35: Two rational tangles, and a pseudoline arrangement. 





Figure 36: Some eight-strand tangles obtained by replacing the intersections of Figure 35 b) 
with rational tangles. 
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12 The Electrical Linear Group Perspective 

In Is], Lam and Pylyavskyy define a Lie algebra zl2n whose associated Lie group EL2n 
is called the electrical linear group. This group has generators Ui (a) for 1 < i <2n and 
a G M, and these generators happen to satisfy the equations 

Ui{a)ui{b) = Ui{a + h) (5) 

Ui{a)uj{b) = Uj{b)ui{a) (6) 

when \i — j\ > 1, and 

Ui{a)uj{b)ui{c) = Uj{bc/{a + c + abc))ui{a + c + abc)uj{ac/{a + c + abc)) (7) 

when |i — j'l = 1 and the denominators make sense. This group is isomorphic to the 
symplectic group Sp2n- One isomorphism goes as follows: let Cjj denote the n x n 
matrix with a 1 in the ith. row and jth column, and zeros elsewhere. Then in block 
notation, we will take 



U2i^i{a) = L j''j £ Sp2 



J2n 
\ U II 

for 1 < i < n, 

U2i{a) = ( / , \ r ) ^ ^'P2n 

\—a[eii + Cj+i^i+i — ej^j+i — Cj+i^jj IJ 
for 1 < i < n, and 

U2n{a) = ( A e Sp2n- 

\-aenn I J 

Within the electrical linear group EL2n we have the electrically nonnegative part 
[EL2n)>o which is the submonoid generated by the Ui{a) for a > 0. The following 
fact is essentially Proposition 2.4 in fsl. 

Proposition 12.1. For i = ii . . .ii a word in the alphabet {1, . . . , 2n}, let ipi : M-^q — t- 
EL2n be the map given by 

ipi{ai,...,ai) = Ui^{ai) ■ ■ ■ Ui^{an) G EL2n- 

Then the following are true: 

• -01 is injective if and only if Sjj • • • Sj„ is a reduced word in the symmetric group 
S2n+i, where Si denotes the transposition of i and i + 1. 

• When Sij^ ■ ■ ■ si^ is a reduced word for some w G S2n+i, the image of ipi depends 
only on w. 

• Letting E[w) denote this image, {EL2n)>o can be written as a disjoint union 
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This proposition is of interest for a couple of reasons. First, it is a disguised restate- 
ment of the main results of ll] and J2|r| Second, it is completely analogous to similar 
facts about factorization schemes for nonnegative upper triangular unipotent matrices, 
which we now recall. Let U be the group of upper triangular unipotent matrices in 
SL2n+i, and let U>o be the submonoid of totally nonnegative ones. (A matrix is totally 
nonnegative if every minor is nonnegative. We refer the reader to [3 or [9] for more 
information on totally nonnegative matrices, and the generalizations of this notion to 
other semisimple Lie groups.) The group U is generated by elementary matrices of the 
form Xi{a) = J + aej^j+i, where o G M and 1 < i < n. These generators satisfy relations 
analogous to ([HIIT]), namely 

Xi{a)xi{b) = Xi{a + b) 

Xi{a)xj{b) = Xj{b)xi{a) 

when \i — j\ > 1, and 

Xi{a)xj{b)xi{c) = Xj{bc/{a + c))xi{a + c)xj{ac/{a + c)) 

when |i — j'l = 1 and the denominators make sense. The submonoid U>q turns out to 
be generated by Xj(a) for a > 0, just as {EL2n)>o is generated by the Ui{a) for a > 0. 
Moreover, the following analogue of Proposition |12.1| holds: 

Proposition 12.2. For i = ii . . .ii a word in the alphabet {!,..., 2n}, let tpi : M.yQ — )■ U 
be the map given by 

V'i(ai,...,ai) =Xii(ai)---Xi„(a„) G EL2n- 

Then the following are true: 

• -01 is injective if and only if Sjj • • • Si^ is a reduced word in the symmetric group 
5*271+1, where Si denotes the transposition of i and i + 1. 

• When Si-^ ■ ■ ■ Si^ is a reduced word for some w £ S2n+i, the image of ipi depends 
only on w. 

• Letting E{w) denote this image, U>o can be written as a disjoint union U>o = 
This is essentially Proposition 2.7 of |9|. 



Lusztig's proof of Proposition |12.2 also establishes the following strengthening of 



the first bullet point of Proposition 12.2 



""One thinks of the U2i-i (a) as the operations of adding a boundary spike of resistance a at the ith bound- 
ary vertex, and U2i{a) as the operation of adding a boundary-to-boundary edge of conductance a between 
the ith and i + 1st boundary vertices. Relation ([7| corresponds to the Y-A transformation. Injectivity of ijji 
becomes equivalent to the weak recoverability of a certain graph. The first bullet point of the proposition 
is tantamount to the fact that a graph is weakly recoverable if and only if its medial graph is critical. The 
second bullet point follows easily from the relations p][7|) , together with a well-known theorem about reduced 



words in Coxeter groups. The third bullet point has to do with Lemma 11.2 
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Proposition 12.3. If i = ii . . . ii is a reduced word, then the map MLq — )• U 

(ai,...,a/) 1-^ Xii(ai) •••Xi,(a/) 
is injective. 

The image of this map may depend on the choice of the reduced word, however. 
The analogous statement in the setting of electrical networks is true, though it doesn't 
seem to be noted in |8): 

Proposition 12.4. If i = ii . . . ii is a reduced word, then the map MLq — )• EL2n+i 

(ai,...,a/) 1-^ Ui-^{ai) ■ ■ ■ Ui^{ai) 
is an injection. 



This is more or less a disguised version of Theorem 6.7 in the case of linear but 
potentially negative conductances. 

More generally, we can disguise Theorem |6.7| in the following fashion. Let S denote 
the set of all bijections / : M — )• M such that /(O) = 0. For every / G S and 1 < i < 2n, 
let Ui{f) be the bijection M^" — > M^" defined as follows: 

U2i-l{f){vi, ■ ■ .,Vn,Ci,. ..,Cn) = {vi, . . .,Vi_i,Vi - /(Cj), Wj+i, . . . ,U„,Ci,. . . , Cn) 

for 1 < i < n, 

U2i{f){vi,. ..,Vn,Ci,. ..,Cn) 

= (vi,... ,Vn,Ci,. . . , Ci-2, Ci - f{Vi - Vi+i),Ci+i + f{vi - Wj+l), Cj+2, ■ ■ ■ , Cn) 

for 1 < i < n, and 

U2n{f){vi,.- ■ ,Vn,Ci,.. . , C„) = {vi, . . . ,Vn,Ci,.. . ,C„_i,C„ - f{Vn))- 

Note that if / is multiplication by a for some a G M^Q; then Ui{f) is a linear function 
for every z, and the corresponding matrix is the same Ui{a) as above. In the case where 
/ is smooth, Ui{f) will always be a symplectomorphism. We do not assume that / is 
continuous, however. 

Then Theorem 6.7 implies the following generalization of Proposition 12. 4[ 



Proposition 12.5. If i = ii . . . ii is a reduced word, then the map from S' to bijections 
on M?^, sending 

{hi ■■■Ji)^ UiiUi) °---° Uijifi), 
is an injection. 

Interestingly, there is also an analog of this for the setting of upper triangular 
unipotent matrices. For 1 < i < 2n, let Xi{f) be the bijection M^""^^ — )■ M^""*"^ defined 
as 

Xi{f){zi, ...,Zn) = izi,..., Zi-i,Zi + f{zi+i),Zi+i, ...,Zn). 

In the case where / is multiplication by some a G M^Qj the map Xi{f) : M^*^"^^ — )• R^"+^ 
is linear and the corresponding matrix is Xi{a) as above. 
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Proposition 12.6. Ifi = ii...ii is a reduced word, then the map from S' to bijections 
on M^"+-^, sending 

{fi,---,fl) ^ Xi^ifl) ° ■ ■ ■ ° Xi^ifl) 
is an injection. 



This can be proven in the same way as Theorem 6.7, though a bit of ingenuity is 



needed to make the arguments of §9^ and ^ 10 carry over to this context. There are 



probably simpler proofs of Propositions 12.4 through 12.6 



The factorization schemes for totaUy nonnegative upper triangular unipotent matri- 
ces are part of a broader story of factorization schemes for totally nonnegative matrices. 
In light of this, we make the following conjecture, which would generalize results of 
Lusztig |9|: 

Conjecture 12.7. For 1 < i < 2n, and / G S, let yi{f) and hi{f) be the bijections 

yi{f){zi, ...,Zn) = {zi,...,Zi, Zi+i + f{zi), Zi+2, ■■■,Zn) 
hi{f){zi, ...,Zn) = {zi,..., Zi-i,f{zi),Zi+i, ...,Zn). 

If ii . . .ii is a reduced word and so is ji . . -jm, then the map on g'+(2n+i)+m, gg^^^^y^^ 
(/l,...,/,,5l,...,52n+l,A:l,...,A:„)GS'+(2n+l)+„^ 

to 

XiAfl) ° ■ ■ ■ ° Xi^ifl) o hi{gi) o ••• o h2n+l{g2n+l) o Vjiikl) o ' ' ' ° Vimikm) 

should be injective. 

13 Closing Remarks 

We have shown that if F is a circular planar graph, then T is strongly recoverable if and 
only if its medial graph is critical. This extends results of flj and [2], who showed the 
same for weak recoverability. An interesting corollary is that F is weakly recoverable 
if and only if it is strongly recoverable. An obvious conjecture is that this holds for 
general F, not necessarily circular planar. However, it may be necessary to put further 
restrictions on conductance functions (like monotonicity or smoothness). For example, 
monotonicity and smoothness are almost certainly sufficient to guarantee non-linear 



recoverability in the highly non-planar 3-dimensional lattice graphs considered in 10 
by a simple argument unrelated to the one used in the present paper. 

The recovery algorithm presented above makes surprisingly few assumptions on 
conductance functions. This suggests that there might be interesting applications of 
our method. Letting currents and voltages live in finite fields might have interesting 
combinatorial applications. In the linear case, we have shown that the rational map 
from conductances to response matrices remains injective and somewhat well-defined 
when extended from (M>o)" to (C \ {0})". This is a purely algebraic statement, which 
might have some uses. 

We also note that Ian Zemke has generalized some of the results and approaches of 



this paper to infinite networks in 13 . 
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